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MODELING  RADICAL  YIELDS  IN  ORIENTED  DNA  EXPOSED 
TO  HIGH-LET  RADIATION 

J.  H.  Miller,1  W.  E.  Wilson,1  C.  E.  Swenberg,2  L.  S.  Myers  Jr5  and  D.  C.  Charlton4 
'Pacific  Northwest  Laboratory,  Richland,  Wash.,  U.S.A. 

•’Armed  Forces  Radiobiology  Research  Institute,  Bethesda,  Md,  U.S.A. 

Uniformed  Services  University  of  the  Health  Sciences,  Bethesda,  Md,  U.S.A. 

‘Concordia  University,  Montreal,  Quebec,  Canada 

( Received  16  June  1987) 

Abstract — Monte  Carlo  simulation  of  energy  absorption  in  oriented  fibers  of  DNA  is  used  to  model  the 
dependence  of  free  radical  yields  on  the  orientation  of  the  fibers  relative  to  a  flux  of  ionizing  radiation. 

We  assume  a  large  asymmetry  in  the  thermal  conductivity  of  the  fibers  that  permits  rapid  transport  of 
vibrational  energy  along  a  DNA  molecule,  but  not  between  different  molecules.  Based  on  this  assumption, 
our  model  predicts  that  thymine  radical  anions  have  a  significantly  greater  probability  of  undergoing  a 
secondary  protonation  reaction  if  they  are  produced  by  a  flux  of  high-energy  protons  that  is  incident 
parallel  to  the  helical  axis  of  DNA  than  if  they  are  generated  by  a  flux  that  is  incident  perpendicular  to 
the  DNA  molecules.  These  results  are  in  qualitative  agreement  with  experimental  data  on  the  yield  of 
y>,6-dihydroJlwmin-5-yl  radicals  when  samples  of  oriented  DNA  were  exposed  at  77  K  to  neutrons.  _  —  /  . 

Kcua;’jMs'.  H ' ^ k  rjtr  tw  +  r a  \ \  .  r j  L •-  (  L- 

TLr\fl'  T«  ‘i'TrorA.hu  H  rv  r  SAJJr**  0,4 k.  J  LA-''  ^  vJhc 

v  introduction  J  0  mine  am  bn  (T  )  and  guanine  cation  were  the  pnn\ 

cipal  radical  species  present  at  77  K.  However,  when 
trstanding  the  mechanisms  by  which  high-energy  the  neutron  flux  was  approximately  parallel  to  the 
;ed  particles  produce  reactive  chemical  species  in  orientation  of  DNA  molecules,  the  total  yield  of 
gically  significant  materials  is  important  for  radicals  was  reduced  and  a  new  ESR  signal  character¬ 
sing  the  health  effects  of  low-level  exposure  to  istic  of  the  5,6-dihydrothymin-5-yl  (TH‘)  radical 
ing  radiation.  This  is  particularly  true  for  radi-  predominated. 

i  with  high  linear-energy-transfer  (LET)  where  The  most  obvious  difference  in  the  patterns  of 
Nations  between  track  structure  and  target  struc-  energy  absorption  in  the  two  irradiation  geometries 
can  have  a  profound  effect  on  the  nature  and  is  that  charged  particles  that  traverse  the  sample 
bution  of  damage.  The  role  of  nonhomogeneous  nearly  parallel  to  the  helical  axis  of  DNA  have  a 
:sses  in  radiation  chemistry  was  recently  re-  greater  probability  of  multiple  energy  transfer  to  the 
:d  by  Magee  and  Chatterjee.1"  Our  research  in  same  molecule  than  do  particles  incident  on  the 

irea<21  has  focused  on  the  use  of  detailed  Monte  sample  perpendicular  to  the  fiber  orientation.  This 

j  simulation  of  the  spatial  distribution  of  energy  difference  in  the  spatial  distribution  of  energy  absorp- 
rbed  from  charged  particles131  as  a  basis  for  tion  should  have  no  effect  on  the  production  of  free 
rs tending  the  subsequent  evolution  of  chem-  radicals  if  energy  and/or  charge  move  between  DNA 
’  active  species.  In  this  paper  we  report  the  use  molecules  in  a  fiber  as  freely  as  they  are  transported 
lis  approach  to  investigate  production  of  free  along  a  single  DNA  chain.  Hence  the  unusual  orien- 

als  in  oriented  DNA  exposed  to  a  flux  of  tation  dependence  of  radical  yields  in  this  system  may 

energy  protons.  have  important  implications  for  intramolecular 

mples  of  oriented  DNA  fibers  prepared  by  the  energy  and  charge  transfer  in  DNA. 
pinning  technique  have  proven  to  be  a  valuable  Arroyo  et  al.{° 1  speculated  that  the  reduction  of 
•imental  tool  for  investigating  the  production  of  radical  yields  in  the  parallel  irradiation  geometry  was 

adicals  in  DNA  by  both  u.v.  light  and  ionizing  due  to  enhanced  recombination  of  precursor  ion 

tion.|4~,)  The  recognition  that  spatially  cor-  pairs.  This  explanation  seems  plausible  in  light  of 

:d  electron  gain  and  loss  centers  may  be  precur-  recent  measurements  of  conductivity  by  van  Lith  et 

of  DNA  double  strand  breaks11’  has  led  to  a/.110’ which  suggest  that  the  high  mobility  of  electrons 
ved  interest  in  this  experimental  system  for  in-  in  the  ice-like  water  layer  around  a  biopolymer  could 
gating  radiation-induced  DNA  damage.  Re-  enhance  intramolecular  charge  transfer.  Arroyo  et  al. 
y,  Arroyo  et  al.m  used  this  type  of  sample  in  attributed  the  appearance  of  TH’  radicals  in  the 
inction  with  electron  spin  resonance  (ESR)  to  parallel  irradiation  geometry  to  intramolecular  trans- 
tigate  the  production  of  radicals  at  77  K  by  fer  of  triplet  excitation  produced  in  charge  re- 
sure  to  neutrons.  When  samples  were  irradiated  combination.  The  long  lifetime  of  triplet  excitons  at 
flux  of  neutrons  incident  perpendicular  to  the  77K,n>  certainly  favors  this  mechanism.  However, 
al  axis  of  DNA,  ESR  spectra  showed  that  thy-  due  to  the  variation  in  triplet  transfer  rates,"23  triplet 


Understanding  the  mechanisms  by  which  high-energy 
charged  particles  produce  reactive  chemical  species  in 
biologically  significant  materials  is  important  for 
assessing  the  health  effects  of  low-level  exposure  to 
ionizing  radiation.  This  is  particularly  true  for  radi¬ 
ation  with  high  linear-energy-transfer  (LET)  where 
correlations  between  track  structure  and  target  struc¬ 
ture  can  have  a  profound  effect  on  the  nature  and 
distribution  of  damage.  The  role  of  nonhomogeneous 
processes  in  radiation  chemistry  was  recently  re¬ 
viewed  by  Magee  and  Chatteijee.113  Our  research  in 
this  area123  has  focused  on  the  use  of  detailed  Monte 
Carlo  simulation  of  the  spatial  distribution  of  energy 
absorbed  from  charged  particles131  as  a  basis  for 
understanding  the  subsequent  evolution  of  chem¬ 
ically  active  species.  In  this  paper  we  report  the  use 
of  this  approach  to  investigate  production  of  free 
radicals  in  oriented  DNA  exposed  to  a  flux  of 
high-energy  protons. 

Samples  of  oriented  DNA  fibers  prepared  by  the 
wet  spinning  technique  have  proven  to  be  a  valuable 
experimental  tool  for  investigating  the  production  of 
free  radicals  in  DNA  by  both  u.v.  light  and  ionizing 
radiation.14'2’  The  recognition  that  spatially  cor¬ 
related  electron  gain  and  loss  centers  may  be  precur¬ 
sors  of  DNA  double  strand  breaks1*’  has  led  to 
renewed  interest  in  this  experimental  system  for  in¬ 
vestigating  radiation-induced  DNA  damage.  Re¬ 
cently,  Arroyo  et  al.m  used  this  type  of  sample  in 
conjunction  with  electron  spin  resonance  (ESR)  to 
investigate  the  production  of  radicals  at  77  K  by 
exposure  to  neutrons.  When  samples  were  irradiated 
by  a  flux  of  neutrons  incident  perpendicular  to  the 
helical  axis  of  DNA,  ESR  spectra  showed  that  thy- 
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excitons  tend  to  be  restricted  to  the  neighborhood  of 
thymine  that  is  bounded  by  the  nearest  guanine  or 
cytosine  on  each  side.  Hence,  in  DNAs  with  roughly 
equal  G-C  and  A-T  content,  the  range  of  triplet 
migration  is  not  expected  to  be  large.  Ion  quenching 
experiments,,5,  suggest  a  range  of  5-30  base  pairs.  The 
short  lifetime  of  singlet  excitons"41  is  probably  the 
major  factor  that  limits  the  range  of  their  migration. 

In  this  paper  we  present  a  model  for  the  orientation 
dependence  of  TH'  yields  that  is  based  on  an  assumed 
asymmetry  in  the  transport  of  vibrational  excitation 
in  the  oriented  DNA  samples.  Asymmetry  of  phonon 
scattering  rates  in  oriented  DNA  have  been  observed 
using  Raman  spectroscopy;"5'161  however,  these  low- 
frequency  modes  should  not  be  very  effective  in 
activating  proton  transfer  to  T-  if  the  activation 
energy  is  0.2-0.7eV  as  suggested  by  the  work  of 
Graslund  et  al.'S)  Regardless  of  the  mechanism,  our 
model  is  macroscopic  in  that  asymmetric  energy 
transfer  is  approximated  by  the  solution  of  thermal 
diffusion  equations  with  a  coefficient  of  diffusion 
along  DNA  that  is  1000  times  greater  than  the 
coefficient  for  diffusion  perpendicular  to  the  helical 
axis.  This  model  is  developed  in  the  next  section  of 
the  paper  and  our  numerical  results  for  the  prod¬ 
uction  of  TH‘  radicals  by  proton  irradiation  at  0.3 
and  1.0  MeV  are  discussed  in  the  third  section. 
Conclusions  and  the  potential  for  improved  theory 
and  experiment  are  discussed  in  the  final  section. 

THEORY 


affinity,  is  beyond  the  scope  of  this  paper.  Rather  we 
simply  assume  that  the  probability  to  form  T-  in  the 
ith  energy  deposition  event  is  proportional  to  the 
amount  of  energy  e(  absorbed  locally  in  that  event. 
Hence,  the  initial  concentration  of  T-  is  approxi¬ 
mated  by  the  equation 

(4) 

where  gT;  is  the  average  yield  of  T-  per  unit  of 
absorbed  energy  and  <5  (jc)  is  the  Dirac  delta  function. 
Using  this  approximation  in  equation  (3)  we  can 
perform  the  space  integration  and  obtain  the  result. 


£th 

St 


i 

I* 


(5) 


for  the  yield  of  TH‘  relative  to  the  yield  of  its 
precursor  T- .  This  result  is  just  a  weighted  average 
of  the  probability  to  convert  T-  to  TH'  over  all  the 
positions  of  energy  deposition  in  a  DNA  molecule 
from  charged  particles  traversing  the  sample  at  a 
specified  angle  relative  to  the  molecular  orientation. 
The  weighting  factors  express  the  relative  probability 
that  T-  will  be  produced  by  a  given  energy  absorp¬ 
tion  event. 

The  rate  of  proton  transfer  used  in  equation  (3)  to 
calculate  P TH-  depends  upon  the  transient  local  tem¬ 
perature  T(x,  i)  through  the  relation 


k{x,  t)~  A  exp [-Q/Tfx,  r)J,  (6) 


As  suggested  by  Graslund  et  al.,i5>  we  assume  that 
the  major  source  of  TH'  radicals  is  the  reaction 

T'  +XH+-.TH'  +  X,  (1) 

where  the  proton  donor,  XH 4 ,  is  probably  water  of 
hydration.  In  an  excess  of  proton  donor,  this  reaction 
will  be  pseudo  first-order  and  the  yield  of  radiation- 
induced  TH'  radicals  per  DNA  molecule  is  given  by 
the  equation 

yth  =  J  ah 3  dxCT ;  (x,  0)P TH  (x ),  (2) 

where 

*nr(*)- J*  M*.  <)expj^-  J*  k{x,  r'jdr  jdr  (3) 

is  the  probability  that  a  thymine  anion  formed  at 
position  x  on  an  oriented  DNA  molecule,  which  we 
approximate  as  a  long  cylinder  of  radius  b  =  10  A, 
will  be  converted  to  TH’  by  protonation.  The  initial 
concentration  of  T-,  the  rate  of  protonation.  and  the 
duration  of  the  radiation-induced  excess  temperature 
are  denoted  by  CT  ■  (x,  0),  k(x,  f),  and  t^. 
respectively. 

Modeling  the  production  of  T-,  which  probably 
involves  the  decay  of  highly  excited  states  through 
autoionization  channels  followed  by  capture  of  low- 
energy  secondary  electrons  in  traps  with  high  electron 


where  A  is  a  constant  of  the  order  of  vibration  rates 
(10,3/sec)  and  Q  is  the  activation  energy.  If  we 
completely  neglect  intermolecular  transfer  of  vi¬ 
brational  energy  and  make  the  additional  assumption 
that  the  heat  capacity  and  thermal  conductivity  are 
independent  of  temperature,  then  T(x,  t)  can  be 
obtained  by  solution  of  the  one-dimensional  thermal 
diffusion  equation  for  a  cluster  of  point  sources."3141 
This  result  can  be  written  in  the  form 


r(x,r)=rs  +  £rt(x,o,  (7) 

» 

where  7",  is  the  ambient  sample  temperature, 

<»> 

with 


A3 


T 


2Dr> 


(9) 


being  the  parallel  diffusion  relaxation  time  and 


r„ 


y/2Pb:ApC 


(10) 


being  the  initial  excess  temperature  at  xt  due  to 
absorption  of  energy  ik .  The  sample  density,  specific 
heat  and  coefficient  of  thermal  diffusion  parallel  to 
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the  DNA  fibers  are  denoted  by  p,  C  and  Du, 
respectively.  Ultimately  C  and  Dv  must  be  consid¬ 
ered  as  adjustable  parameters  since  we  are  using 
macroscopic  material  parameters  to  describe  a  micro¬ 
scopic  process.  In  addition,  the  effective  specific  heat 
C  contains  a  factor  for  the  conversion  of  energy  from 
electronic  to  vibrational  excitation.  As  a  starting 
point  for  numerical  calculations,  we  use  the  specific 
heat  and  thermal  diffusion  coefficient  of  water 
[4  x  107erg/(g*deg)  and  10_3cm2/s,  respectively). 
The  delocalization  of  vibrational  excitation  by  sub¬ 
picosecond  processes  is  included  in  our  model 
through  the  parameter  A,  which  is  the  initial  width  of 
the  Gaussian  distribution  that  approximates  the  ex¬ 
cess  temperature  induced  by  individual  energy  depo¬ 
sition  events.  We  expect  A  to  be  of  the  order  of  the 
base  stacking  distance  which  is  3,4  A  in  the  B 
conformation. 

The  formulation  of  the  model  thus  far  completely 
neglects  transport  of  vibrational  energy  perpendic¬ 
ular  to  the  helical  axis  of  the  DNA  molecule.This  is 
generally  adequate  for  small  isolated  clusters  of  ex¬ 
citation  where  the  radiation-induced  excess  tem¬ 
peratures  decays  on  a  nanosecond  or  smaller  time 
scale.  However,  when  many  interactions  occur  over 
an  appreciable  length  of  DNA  the  time  dependence 
of  the  excess  temperature  may  be  affected  by  thermal- 
diffusion  perpendicular  to  the  fiber  direction.  Since 
we  are  assuming  that  D,t»D1,  it  is  reasonable  to 
approximate  this  leakage  by  considering  the  loss  of 
thermal  energy  from  a  long  cylinder  of  uniform  initial 
temperature.  Including  only  the  first  term  in  the  series 
expansion  that  is  the  solution  of  this  well-known 
heat  flow  problem,  we  obtain  an  approximation 
for  the  transient  local  temperature  that  is  given  by 

T(x,  /)  =  T,  +  exp^ ~  7-J I  Tk(x,  t),  (11) 

where 


*±  = 


b2 

S.8D1 


(12) 


is  the  perpendicular  thermal  diffusion  relaxation 
time.  If  the  asymmetry  in  thermal  diffusion  is  large 
(we  assume  DLID,f  =  10"’)  and  if  the  initial  width 
parameter  does  not  greatly  exceed  the  radius  of  the 
molecule,  then  t1»ttt  and  the  leakage  factor  in 
equation  (II),  exp{— r/rxJ,  serves  only  to  limit  the 
duration  of  the  excess  temperature  when  many  ex¬ 
citations  occur  in  the  same  molecule.  Hence  it  mainly 
effects  the  yields  calculated  in  the  parallel  irradiation 
geometry. 


RESULTS 

Application  of  the  formalism  presented  in  the 
previous  section  requires  information  about  the  spa¬ 
tial  patterns  of  energy  absorption  in  a  DNA  molecule 
of  specified  orientation  relative  to  the  flux  of  ionizing 
radiation.  In  our  model  this  information  consists  of 


the  positions  of  energy  deposition  events  xk  and  the 
amount  of  energy  tk  deposited  locally  in  each  event 
Energy  transported  away  from  xk  by  high-energy 
secondary  electrons  is  not  included  in  tk,  since  the 
detailed  structure  of  these  secondary  tracks  is  in¬ 
cluded  in  our  Monte  Carlo  simulation  of  the  slowing 
down  of  the  primary  ion.  To  approximate  the  spatial 
distribution  of  energy  deposited  in  a  DNA  molecule, 
we  superimpose  a  long  right-circular-cylinder  of 
radius  b  =  10  A  on  a  computer-simulated  segment  of 
the  track  of  a  high  energy  proton  slowing  down  in 
water.  The  computer  simulation  is  carried  out  by  a 
Monte  Carlo  code  developed  by  Wilson  and 
Paretzke.,3)  The  superposition  of  cylinders  and  track 
segments  is  random  in  all  aspects  except  the  angle 
between  the  axis  of  the  cylinder  and  the  velocity  of 
the  proton.  The  collection  of  energy  transfer  points 
in  this  superposition  of  target  and  track  we  will  refer 
to  as  a  “hit”. 

Once  an  energy  deposition  event  within  a  hit  has 
been  randomly  selected  to  be  the  site  of  formation  of 
T-,  the  transient  temperature  at  that  position  is 
calculated  for  equation  (11).  Typical  results  for  hits 
that  involve  small  and  large  amounts  of  energy 
deposited  in  DNA  are  shown  in  Fig.  I.  If 
^0.001  then  the  leakage  term  in  equation 
(11)  does  not  significantly  effect  the  transient  tem¬ 
perature  at  times  less  than  about  0.1  ns.  However, 
transport  of  thermal  energy  perpendicular  to  the 
DNA  fibers  does  influence  the  excess  temperature 
associated  with  large  hits  at  longer  times  and  substan¬ 
tially  reduces  the  probability  that  protonation  will 
occur.  Unless  otherwise  noted,  all  of  the  results 
presented  in  this  paper  are  for  D±/Dv  =  0.001. 

Assuming  an  initial  delocalization  of  vibrational 
excitation  equal  to  3.4  A,  Fig.  2  shows  the  de¬ 
pendence  of  the  yield  of  TH'  per  T-  precursor  on  the 
activation  energy  for  the  proton  transfer  reaction.  In 
these  calculations  we  have  averaged  over  the  distribu¬ 
tion  of  hits  that  occur  under  a  given  irradiation 


Fig.  1.  Transient  radiation-induced  temperature  changes 
for  two  typical  hits  with  different  spatial  and  energy  charac¬ 
teristics.  For  times  <0.1  ns  temperature  profiles  are  essen¬ 
tially  the  same  for  the  two  levels  of  asymmetry  in  thermal 
conductivity. 
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Activation  Energy  (eV) 

Fig.  2.  Yield  of  TH'  per  T-  as  a  function  of  the  activation 
energy  for  proton  transfer. 


condition.  Typically,  a  sample  containing  on  the 
order  of  103  hits  was  examined.  To  obtain  results  for 
comparison  with  experiment  we  must  further  average 
the  yield  of  TH'  over  a  distribution  of  activation 
energies  for  the  protonation  reaction.  If  we  use  the 
distribution  proposed  by  Grasiund  et  al.,,s>  we  obtain 
the  results  shown  in  Fig.  3.  In  this  figure  the  ratio 
g-mlgr-  is  plotted  as  a  function  of  the  initial  width 
parameter  in  units  of  the  base-pair  separation,  3.4  A. 

Generally,  protonation  of  T-  decreases  ex¬ 
ponentially  with  increasing  initial  delocalization  of 
the  vibrational  excitation  of  the  medium  in  an  indi¬ 
vidual  energy  deposition  event.  This  is  expected  from 
equation  (10)  which  shows  that  the  initial  excess 
temperature  associated  with  any  energy  deposition 
event  is  inversely  proportional  to  the  initial  width  A. 
This  parameter  also  enters  the  calculation  through 
equation  (8)  where  the  effect  of  a  larger  A  is  to 
increase  the  duration  of  the  excess  temperature, 
which  favors  protonation.  Apparently  the  former 
effect  predominates.  In  Fig.  3  we  see  that  protonation 
of  T-  in  the  parallel  irradiation  geometry  is  consid¬ 
erably  less  sensitive  to  the  value  of  A  than  it  is  in  the 
perpendicular  case.  In  the  parallel  case  a  hit  is  made 
up  of  many  more  energy  deposition  events  than  in  the 
perpendicular  case.  This  makes  the  effectiveness  of  a 
hit  for  inducing  protonation  of  T -  less  dependent  on 
the  excess  temperature  associated  with  individual 
energy  deposition  events. 

The  effects  on  the  yield  of  TH’  of  a  ±10°  uncer¬ 
tainty  in  the  orientation  of  the  DNA  fibers  relative  to 
the  proton  flux  is  also  shown  in  Fig.  3.  Results  for  the 


Initial  Width  (base  pairs) 

Fig.  3.  Yields  of  TH"  per  T-  averaged  over  the  distribution 
of  activation  energies  given  in  Ref.  5  as  a  function  of  the 
initial  width  parameter  of  the  model  expressed  in  units  of 
the  base-pair  separation  of  DNA  in  the  B  conformation 
(3.4  A).  Solid  symbols  show  the  effect  of  a  ±10°  uncertainty 
in  the  orientation  of  DNA  molecules  relative  to  the  proton 
flux. 


parallel  orientation  are  more  sensitive  to  this  type  of 
uncertainty  than  are  those  obtained  for  protons 
incident  nearly  perpendicular  to  the  DNA  fibers.  This 
is  a  consequence  of  the  fact  that  the  large  TH'  yields 
in  the  parallel  irradiation  geometry  result  from  the 
relatively  few  hits  where  a  DNA  molecule  is  close  to 
the  axis  of  a  proton  track  for  an  appreciable  distance. 
This  aspect  of  our  model  is  further  illustrated  by 
Table  I  which  gives  the  average  energy  deposited  in 
a  DNA  molecule  by  a  hit  and  the  root-mean-square 
separation  of  T-  from  other  energy  deposition  events 
in  a  hit.  The  latter  quantity  is  a  measure  of  the 
distance  over  which  thermal  energy  diffuses  in  the 
process  of  activating  proton  transfer  to  T- .  Note  that 
this  energy  transfer  distance  is  essentially  equal  to  the 
diameter  of  the  target  in  the  perpendicular  case  and 
is  insensitive  to  a  i  10°  uncertainty  in  this  orien¬ 
tation.  However,  this  level  of  uncertainty  in  orien¬ 
tation  of  DNA  relative  to  a  proton  flux  that  is 
parallel  to  the  fibers  causes  a  large  reduction  in  the 
average  energy  transfer  distance  and  a  correspond¬ 
ingly  large  reduction  in  the  protonation  of  T- . 


Table  I  Characteristics  of  the  interaction  of  oriented  DNA  with  I  MeV 


protons 


Orientation* 

Energy  deposited  (eV) 

Transfer  distance  (A)b 

0 

293 

2240 

0±  10 

150 

474 

90 

62 

21 

90  ±10 

62 

22 

■  Angle  in  degrees  between  helical  axis  of  DNA  and  proton  velocity. 

*  Root-mean-square  separation  between  sires  of  T  and  other  excitations. 


Modeling  radical  yields  in  oriented  DNA 
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The  characteristic  of  our  model  that  is  most  readily 
testable  is  its  dependence  on  the  LET  of  the  radiation. 
This  aspect  of  our  results  is  better  illustrated  by  Fig. 
4  where  the  ratio  of  yields  of  TH'  per  T‘-  precursor 
calculated  with  the  same  model  parameters  in  the 
parallel  and  perpendicular  case,  is  plotted  as  a  func¬ 
tion  of  the  initial  width  parameter  in  units  of  the 
base-pair  separation.  This  ratio  is  an  indicator  of  the 
enhancement  of  the  protonation  reaction  in  samples 
exposed  to  a  flux  of  high-LET  radiation  parallel  to 
the  orientation  of  the  DNA  fibers.  In  this  model,  the 
degree  to  which  protonation  of  T-  is  enhanced  in  the 
parallel  irradiation  geometry  is  strongly  dependent 
upon  the  initial  spatial  delocalization  of  vibrational 
excitation  in  the  DNA.  As  was  discussed  above, 
this  results  from  the  fact  that  protonation  of  T-  in 
the  parallel  case  is  less  dependent  on  high  local 
excess  temperature.  However,  the  increase  in  the 
enhancement  factor  with  increasing  LET  is  essentially 
independent  of  this  initial  width  parameter.  Hence 
experimental  results  at  one  proton  energy  can  be  used 
to  determine  A  for  comparison  between  theory  and 
experimental  at  other  proton  energies.  This  approach 
to  experimentally  testing  the  model  assumes  that  A  is 
independent  of  LET. 

CONCLUSIONS 

Our  model  calculations  show  that  intramolecular 
energy  transfer  following  radiation  exposure  should 
enhance  protonation  of  T-  and  that  this  is  a  reason¬ 
able  explanation  for  the  orientation  dependence  of 
TH‘  yields  that  is  observed  when  a  flux  of  high-LET 
radiation  is  absorbed  in  samples  of  oriented  DNA  at 
77  K.(”  The  amount  of  enhancement  predicted  by  the 
model  is  sensitive  to  the  uncertainty  in  orientation 
and  to  the  stopping  power  of  the  incident  radiation. 
Both  of  these  results  point  out  the  need  to  repeat  the 
expetiments  of  Arroyo  el  al  m  with  proton  irradiation 


Fig.  4.  Enhancement  of  protonation  of  T-  formed  by 
irradiation  with  a  proton  flux  that  is  parallel  to  the  orien¬ 
tation  of  DNA  relative  to  the  yield  of  this  secondary 
reaction  in  the  perpendicular  irradiation  geometry.  Solid 
symbols  show  the  effect  of  a  ±  10°  uncertainty  in  orientation 
of  DNA  molecules  relative  to  the  proton  flux. 


from  an  accelerator  rather  than  a  neutron  source. 
This  will  reduce  the  uncertainty  in  orientation  of  the 
proton  flux  relative  to  the  DNA  fiber  direction  and 
allow  the  enhancement  of  TH‘  radicals  in  the  parallel 
irradiation  geometry  to  be  measured  as  a  function  of 
the  LET  of  the  incident  radiation.  Experiments  of 
this  type  are  in  progress.  Improvements  in  the  model 
are  being  focused  in  two  areas:  (1)  a  more  micro¬ 
scopic  description  of  TH‘  production  with  parameters 
that  are  more  directly  related  to  molecular  properties 
of  the  system;  and  (2)  a  more  complete  theory  that 
includes  recombination  which  will  allow  us  to  predict 
both  the  enhancement  of  TH‘  and  the  reduction  of 
total  radical  yields  in  the  parallel  irradiation  geome¬ 
try.  By  this  combined  experimental  and  theoretical 
effort  we  expect  that  oriented  DNA  will  be  an  even 
more  useful  tool  for  investigating  the  interplay  of 
target  structure  and  track  structure  in  determining 
the  production  of  reactive  species  in  biological  sys¬ 
tems  exposed  to  high-LET  radiation. 
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